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We present the microscopic kinetic theory of a homogeneous dilute Bose con- 
densed gas in the generalized random phase approximation (GRPA), which 
satisfies the following requirements: 1) the mass, momentum and energy 
conservation laws; 2) the H-theorem; 3) the superfluidity property and 4) the 
recovery of the Bogoliubov theory at zero temperature . In this approach, the 
condensate influences the binary collisional process between the two normal 
atoms, in the sense that their interaction force results from the mediation of 
a Bogoliubov collective excitation traveling throughout the condensate. Fur- 
thermore, as long as the Bose gas is stable, no collision happens between 
condensed and normal atoms. In this paper, we show how the kinetic theory 
in the GRPA allows to calculate the dynamic structure factor at finite tem- 
perature and when the normal and superfluid are in a relative motion. The 
obtained spectrum for this factor provides a prediction which, compared to 
the experimental results, allows to validate the GRPA. 

PACS numbers: 03. 75. Hh, 03.75.Kk, 05.30.-d 



1. INTRODUCTION 

The kinetic theory of an interacting diluted Bose condensed gas remains 
still a challenging problem. The situation is indeed not entirely clear even 
though some microscopic description of the condensate dynamics have been 
attempted by many groups EEEI Some of the difficulties are really related 
to the basics of kinetics that Boltzmann developed more than one century 
ago, in order to derive an equation for the particle position and momentum 
distribution in a classical gas. This equation obeys the conservation laws for 
the mass, the momentum and the total energy density as well as the famous 
H-theorem which allows to recover the second principle of thermodynamics. 
Strangely, the extension of that approach to a quantum Bose condensed gas 
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Fig. 1. Feymann diagram of the interaction of two normal atoms (full line) 
mediated by the Bogoliubov collective excitation (wavy line). 



compatible with these conservation laws and the H-theorem has not seemed 
to attract a lot of attention. In addition, the kinetic equation must describe 
the superfluidity, in particular, the persistence of a non zero relative veloc- 
ity between the normal and superfluid when it is below the critical velocity. 
Also, we expect to recover at zero temperature the Bogoliubov theory whose 
validity has been confirmed experimentally -. At first glance, all these re- 
quirements seem logical and well defined but an analysis of the literature 
reveals the difficulties to satisfy them in a unified theory. But, a solution to 
this problem has been recently presented with surprising properties^. 

A kinetic equation has been derived by extending the generalized ran- 
dom phase approximation (GRPA), generally used for describing charged 
particles, to the case of a condensed atom gas. For a system of charged 
particles, the resulting equations describe the binary collision between the 
particles as a process of emission and absorption of a plasmon mediating 
the interaction. This plasmon corresponds to the longitudinal component 
of an electromagnetic wave in addition to the transverse ones. For neutral 
bosons, however, the condensate itself helps the thermal atoms to collide 
by transferring momentum and energy. The vehicle for this transmission is 
the Bogoliubov collective excitation energy with its dispersion relation char- 
acterizing a Goldstone boson (see Fig.l). Surprisingly also, provided that 
the condensate is stable, particle exchanges between the condensate and the 
thermal cloud are excluded since no collisions allow such a process. This 
point can be understood qualitatively by analogy with salt in the water. 
The interaction force between ions of salt is strongly attenuated due to the 
dielectric power of water to generate a counter polarized force. In the same 
way, due to its powerful "dielectric" behavior, the condensate can deform 
its wave function in a coherent way in order to diminish its potential energy 
and attenuate locally the potential interaction force displayed by any excited 
atom. In the GRPA, this effect is so strong that any atom in the conden- 
sate does not feel the presence of an excited atom anymore leading to the 
absence of any collision process. As a result, the particles of the normal and 
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superfluid become invisible to each other and feel only the mean field force 
of the other fluid. This different picture allows to understand more simply 
why a relative velocity between the two fluids persists forever and leads to 
a metastable state characteristic of the phenomenon of superfluidity. This 
metastable state persists in time provided that the relative velocity does not 
exceed the sound velocity. Otherwise, the Landau criterion is not satisfied 
and the gas becomes unstable leading to some particle exchanges again. 

The purpose of this paper is to show how to calculate the dynamic 
structure factor in the GRPA at finite temperature and in presence of a rel- 
ative velocity. Since this factor can be measured in experiment, comparisons 
can serve to validate the GRPA approach. After recalling the basics of the 
GRPA method for the kinetics of a Bose gas in section 2, we show how to 
calculate this factor in section 3. 



2. THE KINETIC EQUATION IN A HOMOGENEOUS GAS 

For any momentum q, we define the density operator p q = X^kPk,q 
where /5k,q = c^Ck+q and Ck (cjj is the destruction (creation) operator. 
The GRPA assumes that the diagonal contribution is the dominant one 
i.e. (/5k, q) = <5q,o n k- For q = 0, the time evolution in a volume V of the 
population average of the atom with mass m in the mode k is given by^ 

^ nk = 12 mV^ Pk ^' ~ ^ k -q'.q')^q' + Pq>(Pk,l' ~ Pk-q',q')) (1) 

q' 

where a is the scattering length. For q 7^ in the GRPA, we neglect con- 
tribution that are quadratic in the non-diagonal operator in order to get an 
approximated linear equation: 
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y XT [™ k ( 2 ~ ^ k > k ' ~~ <*k-q,k') ~~ n k+q(2 - <5k,k' - <Wq,k')]/°k',q (2) 
171 k' 

where ek = k 2 /2m. This equation possesses two types of eigenvalues^ 1) 
the scattering solution with the spectrum given by the kinetic energy trans- 
fer: uj = ek+q — £k! 2) for excitations involving the macroscopic condensate 
mode nk s i.e. for k = k s ,k s — q, the presence of interaction transforms the 
scattering solutions into collective solutions u> = lo q — «7 q . These are the 
zeroes of the dynamic dielectric function (DDF): 

87rcro ks q 2 

K(q, W ) = Mq, ») - ( „ + i0 h^r ,«■),+ «-*<■ (3) 
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where e q = -^/c 2 q 2 + (q 2 /2m) 2 is the Bogoliubov excitation energy, c = 
y/kitani^ /m 2 V is the sound velocity and 

VI \ 1 8vr ° V^ n k _n k+q m 

/C n (q,cj) = 1 > ; — - 5- (4) 

k ' m 2m 

is the DDF of the normal fluid where n^ = (1 — <5k,k s ) n k- The real part and 
imaginary part correspond to the Bogoliubov collective excitation and the 
Landau damping respectively. For a temperature close to zero, there are two 
solutions 1 : o;± ~ ^ ± ef, 7 ± ~ ±Im/C n (q,u;±) C 2 q 2 /e^. 

The substitution of the solution of Eq.(J2J) into (^0) allows finally to get 
the GRPA kinetic equation for a Bose gas™ As explained in the introduction 
we find that dn^/dt = i.e. the condensed particle number is a constant 
of motion and thus no particle is exchanged with the normal fluid. For the 
normal fluid, the populations in each mode obey the balance equations: 



q,k' 



in V 



71"£(ek+q + £k'-q ~ Ck ~ Ck') 



£(q,£k+ q -ek) 

Kc+q n k'-q( n k + !)«' + 1) - n k n k'( n k+q + ^Kc'-q + X )] ( 5 ) 

These show that two normal particles collide through emission and absorp- 
tion of a Bogoliubov collective excitation, called "condenson" (see Fig.l). 
The dispersion relation for an excitation of momentum q is given by the 
pole of the DDF which plays the role of a propagator in ©. For a mo- 
mentum q, this excitation has the energy w q and a decay rate, due to its 
absorption by a thermal atom, given by the Landau damping 7 q . Similarly 
to the plasma case^, these equations are valid provided that the collective 
oscillations are damped i.e. 7 q > 0. Otherwise, they are amplified exponen- 
tially and the Bose gas becomes unstable. The case with instability leads to 
different kinetic equations and is not treated here. But the instable regime 
allows particles exchange between the two fluids. The equilibrium solution 
of © for a normal fluid at rest is given by n£ q = 1/exp [/3(ek — /i)] — 1. In 
this case when the temperature is close to zero, we find that: 

ImK „ (q ,. ) = wln ^___ I _^_ E _j (6) 

which is a positive function for u> = u^ > and negative otherwise. Thus, 
we deduce that the stability condition 7^ > is equivalent to the Landau 
criterion e q > k s .q/m i.e. |k s /m| < c. In other words, the condensate is 
stable if its velocity relative to the normal fluid is smaller than the sound 
velocity. 
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3. THE DYNAMIC STRUCTURE FACTOR 

The dynamic structure factor is defined as: 

{' f'OO J / 

5(q,w)= lim dr e -(q™') ( (r) p _ q ( r + r ')) (7) 

This function is calculated within the GRPA and within the assumptions 
that the normal fluid is in thermal equilibrium at rest. We start with the 
initial condition that for r = the correlation function contains only non 
interacting contribution^! (/9k',- q Pk,q)|r=o = ( n k + l)"-k'<^k'-k,q(l - <5 q ,o) + 
n k n k"5q,o- Through Eq.(5), the evolution of the system will create for r — > oo 
some correlation due to the presence of the interaction between particles. 
Similarly to the plasma case, the resulting dynamic structure factor is pro- 
portional to the inverse of the DDF : 

1 1 mV / 1 \ 

2n exp(— pu!) — 1 Aira \K,(q,u))J 

The relevant contributions are near the poles of /C(q, uj). At temperature 
close to zero, we obtain the approximated expression: 

^' 1 +^^) (9) 



^(q> W ) e q \LU-U^ + i^ UJ - W q + ?7 q 

Eq.© and @ are the main results of this paper and generalize previous 
results "El for any temperature in the presence of a relative velocity between 
the two fluids. Due to the asymmetric dispersion relation, this factor has 
two asymmetric peaks with usual intensity given by <S*(q) = q 2 /2meq . Also 
the presence of temperature gives rise to the broadening of the peaks with 
a width corresponding to the Landau damping. When the relative velocity 
approaches the sound velocity, the asymmetric Bogoliubov frequency ap- 
proaches zero for some non zero momentum for which the Landau damping 
gets closer to zero. Beyond that velocity, in some regions of the momentum, 
both the resonance frequency and the Landau damping become negative. In 
that case, the gas becomes instable leading to non equilibrium dynamics. 
Thus the experimental analysis of the dynamic structure factor provides a 
different way to confirm that the sound velocity is the critical velocity be- 
yond which some instability should appear. In real experiments, the width 
comes mainly from inhomogeneous broadening but, in principle, the asym- 
metry between the peaks must be observed. For this purpose, we need to 
create a situation with a non zero relative velocity^. Knowing the spatial 
density and velocity distribution of the condensate, we can use the local 
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density approximation to calculate the line shape of the Bragg excitations 
and test the validity of the results ™*l 

To conclude, the GRPA formalism confirms previous results established 
for the dynamic structure factor. Moreover, it allows to calculate this quan- 
tity from not only for an equilibrium distribution function but also for a 
non equilibrium one. Of special interest is the metastable state of superflu- 
idity from which an explicit expression has been obtained for temperature 
close to zero and which displays the transition from metastability to insta- 
bility. This transition is dictated by the Landau criterion which delimits the 
region in which the dynamic structure factor is well defined. Beyond that 
region this function is not stationary anymore due to unstable dynamics. 
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